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Invited  Paper 


Theory  of  Second  Harmonic  Generation  in  presence  of 
Diffraction,  Beam  Walk-Off  and  Pump  Depiction 

Shekhar  Guha  and  Lconel  P.  Gonzalez 
Materials  and  Manufacturing  Directorate 
Air  Force  Research  Laboratory 
Wright  Patterson  Air  Force  Base,  Ohio  45433-7702 
email;  ahekhar, guha@wpafb.af.mil 

ABSTRACT 

Integral  expressions  tor  the  pump  and  generated  fields  are  presented  here  for  the  case  of  second  harmonic 
generation  of  a  focused  Gaussian  pump  beam  incident  on  a  nonlinear  crystal  .  The  birefriflgfcnt,  walk-off  of  the 
generated  beam  and  the  effect  of  pump  depletion  are  included  in  the  theory. 


1,  INTRODUCTION 

The  theory  of  second  harmonic  generation  (SHG)  of  a  focused  Gaussian  beam  in  the  presence  of  beam  walk-off 
presented  by  Boyd  and  KleinmantJ-  has  been  widely  used  in  the  nonlinear  optics  community.  A  limitation  of 
their  method  is  the  exclusion  of  pump  depletion  which  makes  it  inapplicable  to  the  case  of  high  conversion 
efficiency,  The  numerical  method  developed  by  Smith  C3-1  and  recent  work  by  Wang  and  Weiner  (3:i  and  by 
Kasamatsu.  Kubomnra  and  Kan  extend  the  theory  of  Boyd  and  Kleinman  to  include  the  effect  of  pump 
depletion,  We  present  here  an  alternate  technique  to  calculate  the  spatial  distribution  of  the  pump  and  the  3HG 
fields  in  presence  of  diffraction,  linear  absorption,  phase  mismatch,  beam  walk-off  and  pump  depletion.  The 
field  distributions  are  expressed  as  multiple  integrals,  which  are  considerably  simplified  for  the  case  of  collimated 
beams,  i.eM  when  the  spread  of  the  pump  beam  due  to  diffraction  inside  the  crystal  can  be  ignored.  This  is  often 
the  case  when  high  peak  power  pulsed  lasers  are  used  as  pump  beams  because  the  pump  is  intentionally  kep; 
collimated  through  the  nonlinear  optical  medium  to  avoid  damage  to  the  material 

In  the  theory  presented  here  the  incident  pump  beam  is  assumed  to  be  focused  at  the  crystal  center  -  this 
restriction  can  he  easily  removed  later.  The  aim  of  this  work  is  to  describe  the  general  approach  taken  and 
to  provide  expressions  for  the  electric  fields  expanded  up  to  fifth  power  in  the  second-order  nonlinear  optical 
coefficient,  d.  The  method  described  here  is  an  extension  of  an  early  work  on  optical  parametric  oscillation 
(OPO)1-'-"-  and  can  be  react ily  applied  to  other  nonlinear  optical  processes  of  ennenl,  interest,  such  as  pump 
resonant  SRO  or  resonant  second  harmonic  generation  (SHG).  As  pointed  out  in  ref-3,  the  Green's  function 
method  is  easily  extended  to  the  time  dependent  case  as  well.  However,  a  limitation  of  the  perturbative  approach 
used  here  is  that  multiple  integrals  need  to  be  evaluated,  which  becomes  computationally  time  consuming  when 
six,  eight  or  higher  dimensionni  integrals  are  involved. 

In  the  following  sect  ions,  the  derivation  of  the  relevant  wave  equations  is  provided  first,  to  set-  up  the  starting 
equations.  Then,  the  Green's  function  method  of  solution  in  terms  of  a  power  series  in  the  nonlinear  coefficient 
is  described,  and  the  expressions  for  the  pump  and  the  SHG  beam  electric  fields  are  evaluated  for  a  few  terms. 


2+  THE  STARTING  EQUATIONS: 

The  notations  for  electric  fields  and  nonlinear  polarizations  used  here  will  generally  follow  the  convention  adopted 
in  the  Handbook  of  Nonlinear  Optics^  {using  the  Si  system  of  units).  The  subscripts  p  and  s  denote  the  pump 
and  the  SIIG  fields,  respectively.  The  fields  are  assumed  to  be  linearly  polarized,  and  only  the  scalar  components 
of  the  electric  fields  are  considered  in  the  nonlinear  interaction.  The  equation  16,4.0  in  Quantfum  Electronic#*'7- 
is  taken  as  the  starting  point  of  the  analysis  lie  re: 


8E  Q* E  d2 

^m+tlo^+tloWp^L 


(i) 


i 


Nonlinear  Frequency  Generation  snd  Conversion  Materials,  Device*,  end  Apclicallons  VI  soiled  by  Peter  E  Powers 
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(2) 

(3) 


The  electric  fields  at  the  pump  and  the  SHG  frequencies  and  w#.  respectively)  are  written  as 

Ew'r{r,  t)  =  Epfrje-^*  +  cr. 

Er-{v,t)=Bf(t)e~iiJ’t +cc 

and  the  total  nonlinear  polarization  as 

Pnl  [r,  t)  =  (P£’t(r)e-^*  +■  cc)  +  (PfaMe-'"-'  +  «)  (4) 

Following  ref.  6  (chapter  2)  -we  assume  that  the  nonlinear  polarization  terms  at  the  pump  and  the  SHG  frequencies 
can  be  written  as 

Ki(r)  =  4^£,(r)E;(r)  (5) 

KlW  =  2de0^(r)a  (6) 

Inserting  equations  213,4>5  and  6  in  equation  1,  and  collecting  terms  oscillating:  at  frequencies  ^p  and  uX7  we 
obtain, 

(V4  +  ^1)£p[r)=-^XiW  (7) 

( Va  +  kl  )E,(r)  =  -ftwXW  (8) 

where 

ftp,  =  A:a  -T-  =  Jca  4-  (9) 

' 

In  the  above.  It  has  been  assumed  that  in  equation  1,  £PiJS  =  Mid  that  the  wave  vectors  are  AyiT  =  nPia^P|i/c 

where  np  and  denote  the  refractive  indices  of  the  medium  at  the  two  frequencies. 


3.  GREEN’S  FUNCTION  SOLUTION 


The  solution  to  the  inhomogeneous  differential  equations  7  and  8  above  is  obtained  using  the  Green's  function 
method: 


&*(*} : 


S± f 
4*  / 


,3>lSpj  p 

&  Rn  KW 

(10) 

*  R. 

(ii) 

where  Rp  and  i?*  denote  the  magnitudes  of  the  displacement  vectors  connecting  the  points  r  and  r  for  the  pump 
and  the  SHG  fields  respectively,  The  pump  beam  is  assumed  to  propagate  in  the  medium  as  air  ordinary  wave, 
so  that 


rp  =  I  r  -  r  '  | _ 

=  \/(jt  -  zH)3  +■  (v  ~  ^)s  +  ”  i')S  {12) 

If  die  SHG  lit  Id  propagates  through  the  nonlinear  crystal  as  an  extraordinary  wave,  say  with  a  walk-off  angle 
denoted  by  p,  and  if  the  Cartesian  coordinates  are  chosen  such  that  the  propagation  vector  of  the  SHG  field  lies 
in  the  x,  z  pl&ne  and  is  at  angle  p  with  the  i  axis,  we  b&ve  W 

R*  =  -  pz)  -  (s'  -  pa')}3  +■  in  -  y')‘l  +  (z  -  s')3  (13) 

For  small  absorption,  equation  11  can  be  re-written  as: 

kPl  =  kp  +  kBl  *k.+i^  (14) 

where  Op;i  =  denote  the  the  linear  absorption  coefficients.  The  Fresnel  approximation 

«*  =  (*-  *')  +  -  *T  +  (»  -  ^'}2}  (15) 
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and 


R‘  =  {z~  2>)  +  2{^7)[{Xa  - ^  +{v~ v')3} 


(16} 


is  used  in  the  exponents,  with  xa  =  x  —  p[z  —  z*}  and  the  approximation  Rv,Rh  =  (z  —  z1)  is  used  in  the 
denominators-  Equations  10  and  11  then  reduce  to 


and 


c,  t  W&jp  f  dx'dtfdz* 

E’[r)  =  -^rj 

eiV*-*')-nf ft-O+j^ W)s  W,  f 

NL  ^  f 

r  f-',  f  dr/dy'dz1 

EAt)  =  -—  j  ~ — 

^ihM{z-z  -  — -  )+  a j a, _*£ , ,}  [ -#  JL  — {i— t' ^.i  j 


{17) 


(18) 


The  incident  pump  beam  is  assumed  to  have  Gaussian  cross-sections  and  to  be  focused  at  the  center  of  a  nonlinear 
medium  of  length  if.  Assuming  the  origin  of  the  coordinate  system  to  be  at  the  incident  face  of  the  medium,  the 
amplitude  Ej,u,(r)  of  the  incident  pump  field  can  be  written  as 


E"tti  { r )  = _ ^£5 _ e*fcM  -  ocp*/2  - {x2  -  y2 )/{ ( i+ai(  a  -  i )  /&* } 


(19) 


where  &p  =  denotes  the  con  focal  parameter.  In  the  first  order  approximation,,  ire.,  when  pump  depletion  is 

ignored,  it  is  assumed  that  the  pump  field  Ep  remains  equal  to  £p°\  a  nonlinear  polarisation  at  frequency  ^  is 
generated  through  equation  6,  which  gives  rises  to  an  electric  field  (dEll]  )  at  frequency  y,  through  equation  11. 
This  generated  SHG  field  interacts  with  die  pump  field  to  generate  a  nonlinear  polarisation  through  equation 
o.  which  generates  (through  equation  10}  a  component  of  the  electric  field  at  the  pump  frequency  (denoted  by, 
say,  Following  this  procedure,  a  series  of  SHG  and  pump  field  terms,  respectively  odd  and  even  powered 

in  d ,  can  be  obtained.  That  is,  in  general,  the  SHG  and  the  pump  fields  can  be  written  as 


EI  =  dE}1!  4-  dJEj3>  +  dPE® > . . 

(20) 

Ep  ~  E®*  +  d?EW  +  +  <fE^ . . . 

(21) 

Substituting  these  scries  in 

the  equations  a  and  6,  the  nonlinear  polarization  terms  can 

also  be  expanded  in 

power  series  In  d: 

T*>,  -  .  p^s)  .  W*' 

n VL-  -  FNL  +  FSL  FNL  ■  ■ 

(22) 

_  pe,:3>  +  pp^  ,  -pfm 
hnl  ~  *nl  +  fnl  +  fnl  -  ■  ■ 

(23) 

where 

(24) 

HI  =  4d3c0£^>^2> 

(25) 

(26) 

and 

(27) 

Ht=^&  ^ 

(28) 

HI  =  *Ptv{Ef>Efr  +  E™E™'  +  E^E?‘) 

(29) 
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Using  equations  19  and  24.  and  substituting  in  IS,  the  finst  term  on  the  right  hand  side  nf  eqn-  20  can  be 
determined.  Using  this  expression  for  dE;1j  and  eqn.  10  in  eqn.  27.  the  {FE^  term  in  eqn.  21  can  be  evaluated 
from  eqn.  17.  With  this  tfPE'p' 1  and  the  Es  from  eqn.  19  substituted  in  eqn.  25,  equation  17  can  be  used  for  the 
evaluation  of  the  tflE';*-  term  in  Eqn.  20,  and  so  on.  Before  embarking  on  the  explicit  determination  of  each  of 
these  terms,  in  the  next  section  several  parameters  axe  defined,  and  a  new  notation  is  introduced  which  allows 
all  the  apatial  integrations  to  be  performed  in  dimensionless  units. 

3*1*  Notation 

The  parameter  =  l/hp  where  t  denotes  the  length  of  the  nonlinear  medium  is  used  to  characterize  the  amount 
of  focusing  of  the  pump  beam.  The  phase  mismatch  between  the  nonlinear  polarization  beam  and  the  generated 
SHG  wave  is  denoted  by 

=  2k$,  —  k #.  (30) 

Similarly,  the  mismatch  in  linear  absorption  between  the  nonlinear  polarization  and  the  SHG  beam  is  proportional 


Aq  =  otp-  a,/ 2. 

In  addition,  the  following  parameters  are  defined: 

(31) 

p  =  x/wop  q  =  y/wap  t  =  zfl 

(32) 

"hi 

Ml 

i 

III 

& 

i 

III 

iS 

[33) 

pa  =  xJmop  U-q  =  ptfwb' p 

(34) 

where  Wi  is  a  dimensionless  parameter  which  provides  a  measure  of  the  amount  of  walk-off  of  the  SHG  wave  as 
it  propagates  through  the  crystal.  Also,  defining  the  variables  ap  and  dp(t}  as 


%  —  1  — 

=  +  2*&* 


the  incident  pump  field  can  be  expressed  as: 


E^Ur)  =  —^Leikylt 
p  ■  '  di[t) 


(35) 

{36) 

{37} 


4,  THE  EXPRESSIONS  FOR  THE  ELECTRIC  FIELD  TERMS 
4+T  The  first  order  SHG  field  term 

Using  Eqns.  L8_L9  and  24.  the  first  order  5]gnas  field  term  dE^  is  obtained  (after  integrating  over  the  transverse 


coordinates ) : 
where 

and 


dBj^r)  =  KEpafPi  r) 


2ijrdiEp0 


fP( r)=e<“— fdt: 

Jit 


eft  At- Aeytii  e -(P, 1 !+S*)  i^s^TTT 

{^i)<fe{Mi) 


Here  pa  =  p  —  W\  (f  -  ti),  The  function  d2(i,  i,)  is  defined  as 

dafMi)  =  di(ii)  ■+  2i^p(i  -  t]), 

where  fc  =  2 kp/kt. 


m 

m 

(40) 

(41) 
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4h2+  The  second  order  pump  field  term 

Using  Eqns  38,19.  27  aid  17s  the  second  order  pump  field  term  d2E^'  can  be  obtained  (after  integrating  over 
the  transverse  coordinates}: 

(42) 

np 

where 

fjfH  r)  = 

x 


*a  (Misfit)  =  *  {+  f  +  .1 
d*  (Mi  1*2) 

(44) 

ta(*,*tJa)  =2IVi(»i-*»}  .^l)‘ 

(45) 

and 

I 

cj(i,ii(^}_  2lVi  (*i  is)'  £^t  ^ 

(46) 

Other  functions 

used  above  are  defined  as 

^2a(ti  ti)  =  eft  (*1 T  -  2i^(t  -  ti) 

(47) 

tfe(*ista)  —  ^(M  fa)  +  2(^i(ii)fc 

(48) 

dd(U  t j ,  ^}  =  f/i{ti)¥dg(ti  ,  fa)  Hr  2t|p(t  -  ti  }ds(£i .  £3) 

(40) 

/ 


eiikp-*p/2)u  r  ^  t  e(-< A*-®, ,  m id 

h i 


e  --(p3  i  Is  )*a  -3t^a  -S-ca 
^1  (*4}^s(Mii  *2^-3!  U) 


m 


4.-3.  The  third  order  SHG  field  term 

Using  Eqns-  19,  42,  25  and  IS,  the  third  order  SHG  field  term  d:-' E'!;.:! :  c&o  be  obtained  (after  integrating  over 
the  transverse  coordinates); 

^£f(r)=-(^)^^|2/r(r)  (GO) 

nF 

where 


fW(r)  =  e(*.-*./a)«  f*  i 

2o 

X  f  ’  f  1 

JO  /<} 


where 


g—  {Po.  -i-'V  j  rza--i-  2pu  -Vi  +qj 

^3) 


^(Mii^s)  =  2Wi(ta  -  ig) 


deft, 

di(ia)'^i(ti) 


; ^2n ig)s 


(51) 


(52} 

(53) 
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and 

(54} 

with 

M*1  ^  1  *3)  =  ^i(^l  *  *2n  4-  dlftOdafoi  *i) 

(55} 

sk(Ml!  *2i^s)  =  di(ti)d4(tlTta,ta)  +  **fpO  —  h)ds(ti>t2) 

(M) 

and 

£1  (Mis  *2^3)  =  4ifc^p((  -  €  1  )-cii  (ti  )cfi  (i a ) "' s  —  2d2a  (£1 ,.  i2)d%  ( L 1 1  n  tj  ?  £3) 

(57) 

4.4.  The  fourth  order  pump  field  terms 

i’^-big  Eqns-  38,  42,  50  and  19  in  28  and  17.  the  fourth  order  pump  field  term  d4 E i4*  is  obtained  as  a  sum  of  two 

terms  d*E'^  and  rf"3^^1: 

=  2(^12  1  «  I4  r) 

(58) 

^}(r)  =  (^)a  |  N  |4  fipojjfrr) 

ttr 

(59) 

The  expressions  tor  and  are  given  in  the  next  two  subsections. 

4.4,1.  jg> 

r)  =  e(iK-*,/W  f‘ f‘‘ 

JO  Jo 

x  r~  dtsei-^-«'/w>  Jr'*d(4<('"-A«nt. 


Jo  Jo 

^  -  (ph+*s  Ud  3  t  2ptni  -4-04.1 

(60) 

^1  (t4)dg(t,  tlb  t2'  £*) 

where 

yn  it  4  4  4  +  \  —  ^Sl  ^3p  ^4) 

<*,(!,  W.,M- 

(61) 

Mt.  U.h.ts.U)  -  w, 

(62} 

and 

c4t l  Ml  i  i  *3  >  ti)  —  h  1 2 {#3 (tl  ftJ2l  £3 n  £4 )  —  _^4 (£ .  £ i .  £2 1  ^3 s  }  1 

(S3) 

with 

fla(ti,tait3Sti)  —  (ti  —  tsski)  —  2{fa  —  tjiii{ta)*flfi(fa) 

(04) 

_  i*  *  „  ,  \  _  $3n(*l>f2iW*J 

(65) 

„  /*  +  j  \  -  S4n(t,*i>*2i  ^,*0 

(00) 

ii  tj)  —  (ta  “ 

_(*i  -  t2'rdi{t2:t2,t4i}dj;XhJ^U) 

(67} 

:  ■  h.  £3^4)  =  d j  [t2 ,  tss  ti)dfl((iT <3i  tg,  ti) 

(68} 

^4n(tj  M  £21  tai  t*)  —  —  ti)di(£i)¥fl^(£ij  £2+^21  ts)3 

(69) 
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The  functions  d?  and  d$  are  given  by 

<*7(*ii  ^£3^4.)  ^  ^e(*i^2i^3S^)  +  ^i(ii)4fif& {^21^3^4} 


and 


4.4.2. 


■(*) 


^8  (£j  i  !  ^3  ■  £*  }  =  ^1  [t  1 )  dg  ( tl  7  t2  :  £3  1  ti  )  4-  2 ifp  (f  —  1 1 )  cfy(  1 1  h  i-2 !  £3,  £4  ) 


jW(. r)  =  £<**-*./*&  fl  dtie{-i±k-a,/ 2)itl  r  dfie(iA*-AftM(t 

Jo  J® 


i&k—A.a^lt4 


where 


/  1  t /*’  d*4e<- 
./d  Jn 

^  —  (JJ2  — 2;Wjji-|- 

di (t^jd^tJi,  tSl  £3^4) 


a43(th  £1 .  £3l  *3,  £4)  =  11  2-  3l  ^  , 

“lolMli  £2n  ^3^4.) 


^42  h  f  1  n  £2 ,  £3 1  £4)  =  Wi 


h*h'.U} 


and 

with 


^10 (Mi-  £21  £3 >£4.) 

C43{£j£i|£2j£3>t4)  = 

5t(^li  £21  £4)  —  S(£i  —  £2.)^  (^li  £4)  4-  2(^3  —  £4)^21  (ti .  £3)  d^(ti,t2) 

,fn  ,  ,  n_  gjjp(£ 

.^,.  -L:  2.,  3>  -4}  -  ^(tl|ta)^(t1>fa>fJ^0(tjtl|fetf8it4) 

afti  -  *a12  __  3[t3-£4l2rfai(£]^ar 

^2^1  :£s)  cf-3 { 1 1 !.  ^ 3 :  £4.)" 

The  functions  d$  and  dig  are  given  by 

tfy ( 1 1 £-2 %  £4 )  —  2^4 (ti .  £4 ) ¥  +  *£2 (£1  *£2)^3 ( £3,  £4 ) * 

and 

dio{t,  £1,  £2*  h,  /..a)  s  da(£i ,  £a)fij(£i>  £3;  £4)"  +  2i£p(£  —  £1  }du(£i ,  ta, £3*  £4} 

4.5.  Th#  fifth  Older  SHG  field  term 

Using  Eqns  42,  5B;  50  and  19  in  26  and  IS,  the  fifth  order  pump  field  term  ifiE^  is  obtained  ae  t 
terms  d*^?,  d6^  and  d5^1: 


-  (^)aK  I  «  ^  JWJ?M 

flP 

=  4(~}2^  I  K  |4  ^po/isV) 


and 


d&£<f(r}  =  2{“}afc  |  K  *  £U/^(r) 


where  the  expressions  for  the  terms  and  are  given  below 


m 

[71} 

(72) 

m 

(74) 

(73) 

(76) 

(77) 

(78) 

(79) 

m 

sum  of  three 
(SI) 
(82) 

(S3) 
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4.5.1.  fT 


where 


and 


with 


f&>  (r)  =  e***- -a./2)tt  f  ^  x  eCiit- A=)l*l 

Jo 

x  j  1  f  1 

Ju  Jc\ 

x  f  dt f  dt^iAh~Aa^ltli 
Jo  Jn 

£■—  (fa  i  )&fll  -i-2pfl.ifl  :  +C|1 1 

^i(^a)^i(tsMi2(is  (a ,  tgn  Uy 


si  tu-2,a,i)  (fu(tltl(taitBlt4||,)' 


^51  (^^1,^3 1^3  ^4}  =  Wl 


^I2(^-n  £|: *2,  h-iUth)  ' 


^51  *  1  =  *2 , 1 3  f  U)  —  Wj  2  { 


ffs (£  1 1  ^2 :  i  £4  -  £  J  ^ *^120  f  * i  *  1 1  £j i £  3 ;  *■■] ■■  £5 ) 
d  12  (  £i  £t,  i  £2 ,  *3 1 £ 4  >  £fl-  J 


_£b(£i  !  *2,  *3, *4:  £5)) 


$7  (  h  7  h  i  *3 1  £4  J  h)  =  2  ( t2  —  £3.)  di  (£2  )  "*£4  (£l -  £4  ■  is  ) 

+  2{^  —  ^)^i(^)  d.i(ti^3| i3) 

»(ti,fa,  h, U,ti)  -  2(t2  -  t3)  +  2(fJ  -  ts) 


di(£ls£<L>  h) 
2 


fffl(*i^ii*aii4>£s}  =  2{h  -  ta)2  ,  2fj  *'  -I-  2 ((4  -  £&)J  ,  , . 

a4L*ii £2- £4  4, 

The  functions  dn  and  di2  are  given  by 

*£].]  (£  1 1  £2 f  £.1,  £4 )  —  ^(£2,  £4.  £,;■)  -r  ^(£4 j  £5  *  £-2:  *a) 

and 

^1 2  (£?  £1 :  *2 1  *3  *  £4}  — 

+  £kfP(£  "■*i)dii(£i,£21 


4,5,2,  /<|} 


/2J(rj  =  e^-a^a>“  f*dtie(Ukk  'Aa)lu 

Jo 

X  j''  dtie{-i^~^/2)lti  y*' dtaeti^-A»JM» 

./o  A) 

X  f  *  ’i4*-Q*/a)H*  f**  dj.e(iAt-Aft)jf» 

Ju  Jg 

£  -  (Pd  a  -f-?1  Jajg  -2j>^  fri J  -j-^s-a 
x  — - - - _ 

(^)^i(^)-rfia(£i  £21  *ai  £4-  £5.} 


(84) 

(S3) 

m 

(57) 

(58) 

(59) 

(90) 

(91) 

(92) 


(93) 

k 
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where 


and 


^13^1:  ^2-^3  *4=  ^fl) 

dii[tJi:  t4i  tj)1 


^52  (i=  ^  1 1  t2?  h-  t*) 


ji^Oji^aTj^faj^ Uih) 


f-sa(i,  tut2>hi  £4)  s  Wi1  £5)  +04(*irt3i*3>*4i *s) 

+  ^F^-f]H(fjH02)-2ftlf^!i3lt.1.^)31 

^s(^li  ^2>  til  ^4]  tE)d|4'(£,  fl  1  £a)- 

The  functions  and  dm  are  given  by 

^13(^1?  *4s  ^s)  —  di[ti)dj(t^7t^l  t4,  t&)  4-  is(ti^2:  1  *4**5) 


and 


/; 


(5) 

^3 


where 


and 


4  ik^p {t  -  1 1 ) d ia(t I,t2>ht ft ,  tB ) 


/W(r)  =  e(«--^/2)«  fdhe[iAl  **&' 
Jo 

x  f  ^  dtQ^’  i^k~Ctm/m-  I  * 

Jo  jo 

X  ^  *  dtde^“-a'/2)le< 

Jo  J  Q 

g  “(Pn  ^"i"^  )-aB3-+2,Pfl  ifjaa  "T-C53. 


1 2-53  (i-  t  In^S^i,  £4) 


di  (ts  }di  (/  ■> )  di2  ( ^ ,  t  ^ ,  t‘2 1 1$  >.  £4 ,  £5 ) 
 ta>*3!  £4- £5) 


=  W’l 


dis{i:  tii  ^1^5/ 

diitiJgsih.t^Ui  is) 


dm  fat  is  iaJii £4,^} ' 


c&a(£?  £1 :  ^2i  £sj  it )  =  Wi^{g${t\  jfj,  iji  £4,  £5) 

iHpjt  -  „ £5}2 


^lO^li  *2,  i  J|  faj  ts)dlfl(ti  ^1  ^2  1  h>  k:  £j) 

The  functions  rfio  and  dig  are  given  by 

^u(^]  j  ^2>  ^ii  £-i>  £5}  —  oTi {'t  1  )-ofe ( 4e ,  i j,  £4,  £&)  +  diu(£i>  £j:  £.3 ,  ti-  is} 

and 


J 


£1 1 £4*  £5/  — ■  {£i)dio(f i>  ig.  £3 j  £4*  £5} 

—  i&£p(£  -  tj )dj s(ti,  f£.  £3,  £4,  £5) 


(94) 

(9G) 

(96) 

(97) 

(98) 


(90) 

(100) 

(101) 

(102) 

(103) 

(104) 
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5.  CONCLUSIONS 


Expressions  for  the  pump  and  second  harmonic  fields  are  presented  as  series  of  terms  in  increasing  powers  of  the 
nonlinear  optica!  coefficient.  The  effects  of  bean)  diffraction,  linear  absorption  and  beam  walk-off  are  included. 
An  advantage  of  the  method  presented  hero  over  the  widely  used  numerical  method  may  lie  in  the  increase  in 
the  computation  speed,  especially  in  special  cases.  For  example,  in  the  case  of  a  collimated  pump  beam,  i.e,*  for 
ip  ^  In  the  single*  double  and  triple  integrals  in  the  expressions  for  the  fields  reduce  to  analytical  expressions 
and  the  quadruple  and  quintuple  integrals  reduce  to  much  simpler  integrals  which  can  be  rapidly  evaluated  using 
standard  computers.  Also*  the  method  can  be  extended  to  the  ease  of  pump  beams  of  arbitrary  shape. 
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